Exact relativistic viscous fluid solutions in near horizon extremal Kerr background 
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Realistic accretion disk models require a number of ingredients, including viscous fluids, elec- 
tromagnetic fields and general relativistic corrections. Close to the innermost stable circular orbit 
(ISCO) the latter can be appreciable and (quasi-)Newtonian approximations become unreliable. 
This is particularly true for nearly extremal black holes like GRS 1915-1-105, where the ISCO almost 
coincides with the black hole horizon. To describe the physics close to the ISCO adequately in a 
simplified model we approximate the nearly extremal Kerr geometry by the near-horizon extremal 
Kerr geometry and construct in this background relativistic viscous fiuid solutions with electromag- 
netic fields. We discuss some applications of our solutions and possible relations to the Kerr/CFT 
correspondence. 
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Our understanding of black hole accretion rests on analytic 
models. Marek Abramowicz [Id 



I. INTRODUCTION 

John Wheeler coined the term "Black Hole" more than 
four decades ago [2]- At that time Black Holes (BHs) 
were considered as rather esoteric objects of purely the- 
oretical interest and little physical relevance (for a text- 
book on BHs and a history see Ref. Q). The current 
assessment of the role of BHs in physics has changed dra- 
matically, for several independent reasons. BHs are now 
used in various areas of physics, far beyond the original 
realms of application envisaged in the pioneering first 
decade of BH research, including quantum chromody- 
namics and condensed matter physics. We focus here 
on BHs in astrophysics and, to a lesser extent, on BHs in 
the context of the gauge/gravity duality (also known as 
AdS/CFT correspondence). The case for astrophysical 
BHs is nicely summarized in Ref. Q ■ 

By their very nature, BHs can be seen only indirectly, 
e.g. through X-ray spectra of accretion disks surrounding 
the BH. Therefore, most BHs that have been identified 
through astrophysical observations are not isolated ob- 
jects, but have a binary partner — see for instance the 
list in Ref. which contains 20 confirmed BHs. At the 
top of that list is a specific BH, GRS 1915-1-105 in the 
constellation Aquila, which has particularly interesting 
features. 

The mass of GRS 1915+105 is about 18 solar masses, 
M = (18.1±O.4)M0. This experimental result can be es- 
tablished accurately by considering quasi-periodic oscil- 
lations (QPOs) @. (Earlier independent measurements 
led to compatible results, M = (14 ± A)Mq, but with 
higher uncertainty 0, Q). The smn of GRS 1915+105 
was measured more recently [l^. The extraction of 
the spin from the observational data is very involved. 



The method of Ref. [1, [13] employs a spectral analysis 
of the X-ray continuum, using a suitable accretion disk 
model to determine the accretion disk parameters and 
thereby also the BH parameters. The spin, in particular, 
can be determined very accurately due to its profound 
impact on the BH properties, such as the innermost sta- 
ble circular orbit (ISCO). The ISCO is the orbit closest 
to the BH horizon where matter can remain stationary 
on a circular orbit. The faster the BH spins, the closer 
the ISCO is to the BH horizon. In the limit of maximal 
spinning (extremal) BHs the ISCO coincides with the BH 
horizon. This means that accreting matter in a nearly ex- 
tremal BH can move much deeper into the gravity well as 
compared to non-rotating BHs, and this effect leads to a 
hardening of the X-ray spectrum and a higher efficiency 
for conversion of accreted rest mass into radiation. Both 
of these features are observed for GRS 1915+105. It was 
concluded in Ref. [9| that the dimensionless spin param- 
eter is nearly maximal, a, > 0.98, close to the Thorne 
limit a* < 0.998 [TlJ and to the theoretical upper bound 
a, < 1, with a, = 1 corresponding to an extremal Kerr 
BH. Therefore, the BH GRS 1915+105 is a nearly ex- 
tremal Kerr BH. 

In view of the importance of accretion disks for the 
determination of the BH parameters it would be nice 
to have some exact results available for nearly extremal 
Kerr BHs. Since back reactions are negligible (see 
e-g- [H, [il]), a suitable accretion disk model involves 
a fluid on the background of a nearly extremal Kerr ge- 
ometry. This is still a very hard problem, because viscos- 
ity cannot be neglected and electroma gnet ic fields should 
be considered as well, cf. e.g. [11 [11,13. A simplifica- 
tion that is useful in other contexts is the Newtonian 
approximation, which works well if the accretion disk is 
sufficiently far away from the BH horizon [3, [ll] ■ How- 
ever, for nearly extremal BHs accreting matter can come 
very close to the BH horizon and the Newtonian approx- 
imation is not suitable to describe the dynamics close to 
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the ISCO. Thus, if we are interested to describe the dy- 
namics of accreting matter close to the ISCO we have 
to employ General Relativity. The near-extremality of 
GRS 1915+105 may allow for a perturbative expansion 
in the parameter (1 — a*) < 0.02. To leading order in such 
an expansion the dynamics is described by the extremal 
limit a, = 1. We are specifically interested to describe 
the dynamics close to the ISCO, so we can simultane- 
ously impose a second perturbative expansion, namely a 
near horizon approximation. 

Our starting point is therefore a background geome- 
try that is obtained after a double limit: first, assume 
that the BH is extremal and second, assume that we are 
arbitrarily close to the horizon. The geometry obtained 
in this way is known as Near Horizon Extremal Kerr 
(NHEK) and was constructed by Bardeen and Horowitz 
|17| (the possibility of viewing the vicinity of the extremal 
Kerr BH horizon as a spacetime on its own right was al- 
ready suggested much earlier by the findings of Bardeen 
and Wagoner [l3|). 

In this paper we study exact solutions of perfect and 
viscous fluids on the NHEK background, and also include 
a discussion of electromagnetic fields. 

We consider first the standard scenario of a fluid with 
timelike velocity, u°'Ua = —1, which can describe mas- 
sive as well as massless particles, depending on the equa- 
tion of state. We construct exact perfect fluid solutions 
on the NHEK background for a polytropic equation of 
state. However, the Ansatz that u"" be timelike pre- 
supposes that there is a local heat bath, which defines 
a preferred rest frame. Temperature arises here as the 
zero component of a timelike vector, and transformation 
to a moving frame introduces a Lorentz factor. If the 
boost between the heat bath and the reference frame ap- 
proaches the speed of light the Lorentz factor becomes 
singular, and one should not consider temperature as the 
zero component of a timelike vector, but instead consider 
a lightlike vector. The only natural way to do this is by 
demanding u°'Ua — 0. 

The dynamics of a fluid in the NHEK geometry is 
mapped to the dynamics of a fluid at the extremal Kerr 
BH horizon. Particles moving on stationary orbits at 
the BH horizon have to move at the speed of light. The 
same consideration applies to the "heat bath", because 
there is no meaningful way to define a local rest frame at 
the BH horizon. Everything that remains stationary at 
the extremal Kerr BH horizon necessarily moves with the 
speed of light, and the caveat in the previous paragraph 
applies. Thus, we consider a "lightlike heat bath" , in the 
sense that u°'Ua = 0. We call a fluid with the property 
u'^Ua = "null fluid" (not to be confused with a timelike 
fluid, u°'Ua = — 1, that describes lightlike matter, like a 
photon fluid). 

We construct a family of exact solutions for a viscous 
null fluid on the NHEK background in the presence of 
electromagnetic fields. For a given velocity profile we 
predict uniquely (up to an overall rescaling) the viscos- 
ity function r/. The class of velocity profiles that we en- 



counter has very special properties, and we describe them 
in detail. To address stability issues we discuss linearized 
perturbations and establish some rigidity results, which 
point to the linearized stability of our exact solutions. 

Besides the potential phenomenological interest re- 
lated to observations of GRS 1915-1-100, there are various 
purely theoretical motivation to consider the NHEK ge- 
ometry. It is an interesting geometry on its own right and 
displays several remarkable geometric properties, some of 
which we review in this work. Uniqueness and stability 
results were established recently 19, 20, 21]. Moreover, 
the NHEK geometry was exploited in the context of the 
Kerr/CFT correspondence [22|], which predicts that any 
extremal Kerr BH is dual to a certain conformal field 
theory (CFT). This conjecture has engendered a lot of 
recent interest (cf. e.g. H, [H [H [11 [H [H, [H [s^, 
im, m m, m, [M m, llgyrand it couW be rewarding to 
apply our results also in this context. 

This paper is organized as follows. In section |TT] 
we review salient features of the NHEK geometry and 
geodesies in its background. In section IIIII we construct 
timelike perfect fluid solutions for arbitrary polytropic 
equations of state. We mention difficulties with the in- 
clusion of viscosity. In section IIVI we construct null per- 
fect fluid solutions. We are able to lift them to viscous 
null fluid solutions and to include certain electromagnetic 
fields. In section |V] we address perturbations around our 
exact solutions. In section IVTl we provide a discussion of 
our results and put them into the perspective of two per- 
tinent branches of literature: BH astrophysics and the 
Kerr/CFT correspondence. 

Before starting we mention some of our conventions. 
We use signature — , -f . We define the Riemann 
tensor as R%cd a^F^d - ddT'^bc + F^dF^e - F^f.cF'^de 
and the Ricci tensor as Rm :— R'^tad, where we employ 
the usual Einstein summation convention. We choose the 
sign of the e-tensor so that etre4> > 0. Symmetrization 
of indices is defined by V(aUfc) := ^ (VaUfc -t- WbUa)- We 
employ natural units c — G — 1. 



II. NEAR HORIZON EXTREMAL KERR 

In this section we review the features of Kerr and 
NHEK most relevant to our work. 



A. Kerr geometry and ISCO 

The Kerr metric in Boyer-Lindquist coordinates is 



ds' - - 1 



2Mr 



]dt^ 



AMra ain^O 



dt d(f> 



p2 ) 

2Mra^sm''9\.^,. , ,2 , 2 ..2 



-) sin^ ed<l3^ + ^dr^ + p^ dO^ 
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with 
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A :=r2(l 



2M, 



(2) 



It is determined by two physical parameters: the mass 
M and the angular momentum J = aM. Often it is 
convenient to employ the dimensionless Kerr parameter 
a, = a/M instead of the angular momentum. With no 
loss of generality we consider only positive angular mo- 
mentum, J, a, a, > 0. The Kerr BH event horizon 



M(l + ^l-al) 



(3) 



exists only if the inequality < \ holds. If it is satu- 
rated, 



a, = 1 



(4) 



the Kerr BH becomes extremal. If it is violated, a, > 1, 
the Kerr geometry contains a naked singularity. Physical 
processes like accretion of matter cannot convert a Kerr 
BH into a naked singularity if the cosmic censorship con- 
jecture is true |38|. This suggests that Eq. (g]) provides 
an upper bound for the dimensionless Kerr parameter — 
and hence also for the angular momentum in units of the 
BH mass.^ 

To determine the ISCO we consider geodesies of time- 
like test particles in the plane 9 — 7r/2, cf. e.g. [i^]. The 
radial velocity is determined by 



V 



with the effective potential 



cff 



E 



N 



r 



2r2 



(5) 



(6) 



Circularity of the orbits requires V'^'^ = and 
dV^^ /dr = 0. These conditions allow to solve for E and 
L in terms of r. For the ISCO additionally d^V^ /dr"^ = 
must hold because this condition separates stable from 
unstable orbits. This condition allows to solve for the 
radius of the ISCO r = Tisco: 



M 



= 3 + Va?2+3a2 



T \/(3 + a; + 2^x2 + 3a2)(3 - x) (7) 



^ This point of view was recently challenged by Jacobson and 
Sotiriou [39| who found that a BH could over-spin by dropping 
into it a test particle with suitable spin and/or angular momen- 
tum, a procedure which requires a careful finetuning. However, 
as pointed out in their paper, their analysis does neither take 
into account corrections from gravitational radiation nor from 
self-force effects, and such effects could spoil the finetuning re- 
quired for a violation of the cosmic censorship conjecture. 



where the upper (lower) sign refers to co-rotation 
(counter-rotation) and 



l + (l-a2)i/3 [(l + a,)i/3 + (i„a^) 



1/3 



(8) 



For — !■ we get a; = 3 and the Schwarzschild result 
^isco ~ is recovered. In the extremal limit we get 
X =\. The CO- rotating ISCO lies at r'^^^o — ^"^^ thus 
coincides with the BH horizon tqh in ([3|l: 



, (extremal Kerr) = r^/f (extremal Kerr) = M (9) 



B. Near horizon extremal Kerr geometry 

If a* ^ 1 for many applications it is sufficient to ap- 
proximate the Kerr solution by the simpler Schwarzschild 
solution, which is obtained from ([1]) in the limit a, 0. 
Nearly extremal Kerr BHs, a» w 1, do not allow a very ac- 
curate Schwarzschild approximation in the near horizon 
region close to the ISCO. In that region, however, they 
allow for a different approximation in terms of the the 
NHEK geometry introduced by Bardeen and Horowitz 
[l3| . The NHEK geometry is obtained from the Kerr ge- 
ometry (m as follows. One introduces the dimensionless 
coordinates 



2M 



XM 



M 



and takes the limit A ^ keeping i, f, 
yields the NHEK geometry^ 

dp + dp 



b, 9 fixed. This 



ds^ = Af2(l + cos^e) f 

\ r 



4sin^ 



(1 + C0s2 



dt.2 
f 



d9' 



(11) 



The spacetime is no longer asymptotically flat — for in- 
stance, at 6* = the spacetime along the axis is AdS2. 
Thus, one should not think of the NHEK geometry as an 
approximation to Kerr in the same way as one considers 
Schwarzschild as an approximation to Kerr. Rather, the 
geometry ([Til) describes the (infinite) throat geometry of 
an extremal Kerr BH, so anything that happens within 
this geometry is related to physical processes at (or very 
close to) the event horizon. 

The coordinates used in the line element pT|) do not 
cover the full NHEK spacetime. Global coordinates, 
which we shall call again t, r, 9, (p (at the minor risk of 
confusion with the Boycr-Lindquist coordinates), are ob- 



^ In our conventions the relation J = holds. 
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tained through the coordinate transformation 



f ~ {\/ 1 + cost 
i — fyl + r^ sint 



hi 



cos t + r sin t 



1 + VTTr^sint 



(12) 
(13) 
(14) 



The NHEK metric in these gfobal coordinates is then 
given by 



ds"^ = M2(1 + cos^ei) ( - (1 + r^) dt^ + J 



dr^ 



(1 + C0S2 



{d<l) + r dtf + d9^^ (15) 



The range of the coordinates is i £ (— cxd,cxd), r G 
(— cxd,cxd), (p S [0, 27r), 6 E [0,7t]. Note that 9 maintains 
its original role as polar angle since it does not appear 
in any of the coordinate transformations that led from 
Kerr ^ to NHEK The radial coordinate r is not 

restricted to be positive because no singularity is encoun- 
tered at r = 0. 

We conclude this brief review with some geometric fea- 
tures of the NHEK spacetime (|15|) and refer to [T?| for 
a more complete discussion. The t = const, hypersur- 
faces are spacelike globally, so there are no closed timelike 
curves. For sufficiently small radii, < 1/3 the Killing 
vector dr is timelike. It is also timelike if the inequal- 
ity 2sm9 < 1 + cos^d holds, which is saturated for the 
critical polar angle 6'crit. = 0.82 (about 47.1°). If neither 
of these inequalities hold dr may become spacelike, just 
like in the ergoregion of the Kerr BH. The plot Fig. [T] de- 
picts the norm squared of dr- We discuss Killing vectors 
m more detail in subsection [HCl below. While NHEK 
does not exhibit any BH horizon, there is still a sense 
in which one can label a locus in the NHEK geometry 
with the attribute "horizon" , namely the horizon of the 
Poincare-like patch described by the line-element (fTTj) at 
f oo.^ If A is finite in pop the limit f — > cxd indeed 
corresponds to the locus of the horizon. In global coor- 
dinates this horizon is mapped to arbitrary finite values 
of r (if t is infinite) or to r = oo (if t is finite). Thus, in 
global coordinates there is not just a single value of the 
radial coordinate r that would correspond to the original 
BH horizon, but rather it spreads out through the whole 
NHEK geometry. 



^ The local equivalence between the Poincare-like horizon of 
NHEK and the event horizon of extremal Kerr has an analog 
in the simpler case of the extremal Reissner-Nordstrom BH and 
its near horizon limit AdS2 x S'^, see Ref. [4l| for a recent dis- 
cussion. We thank Tom Hartman for discussions on these issues. 




FIG. 1: Norm squared of Killing vector dt. The ("ergo"-) 
region with spacelike dt is plotted in gray, the region with 
timelike dt is plotted in white. The vertical critical line is 
at r = l/\/3- The horizontal critical lines are at 9 ~ 0.82, 
9 ^ 2.32. 



C. Killing vectors, invariants and geodesies 

The Kerr geometry ([T|) is stationary and axially sym- 
metric. Its Killing vectors 



^0 ^ dt fci = 



(16) 



persist in the NHEK geometry p5)) . The latter has an 
enhanced isometry group [13], namely 5'L(2,R) x C/(l), 
whose algebra is generated by the Killing vectors ko , ki 
and two additional Killing vectors k± defined by 



(^rdt-i(l + r^)dr+d^) (17) 



Generic Killing vectors are neither spacelike nor timelike 
globally. The only exception is fci, which is spacelike 
globally. Scalar fields that are invariant with respect to 
all four Killing vectors can only depend on the coordinate 
6. Vector fields that are invariant with respect to all four 
Killing vectors must be of the form 

[At = rA^ie), Ar = 0, = Ae{9), A^ = A^{9)\ (18) 

Since both the Kerr and the NHEK geometry solve the 
Einstein equations all scalar invariants constructed from 
the Ricci tensor vanish. There is an interesting scalar 
invariant that is non-vanishing for the Kerr geometry, 
namely the Chern-Pontryagin density (also known as in- 
stanton density). It is defined by 



CP := 



1 



.cde 



acd 



(19) 
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and for the Kerr geometry it is given by 
96aM^r cos 9 , 



CP = 



(r2 + cos^fl^G 



Sa^ cos26l) (Sr^ - cos^ 



(20) 

The NHEK geometry also has a non-vanishing Chern- 
Pontryagin density 



CP 



96 cos ( 



M4(l + CO& 



(l-3cos26') (S-cos^e) (21) 



Consistently, the NHEK result ([2T|l coincides with the 
Kerr result (|20p in the near horizon extremal limit a — 
r — M. In both cases the Chern-Pontryagin density inte- 
grates to zero so that the instanton number vanishes. In- 
terestingly, the Chern-Pontryagin density ([2T|l of NHEK 
changes sign at 9 = arccos(l/\/3) ~ 0.955, an angle we 
shall encounter again below. As a consistency check we 
have also compared all other polynomial scalar invariants 
of NHEK and Kerr in the extremal limit and found that 
they always coincide. 

We review now geodesies of timelike test particles in 
the plane 9 = t:/2 ^Yc\. The radial velocity is determined 
by ([5|) with the effective potential 



V 



) - 8LEr 



(22) 



We have normalized the 4- velocity to u°'Ua = — 1. The 
constants of motion are given by 



E = u 



2Lr 



We generalize now to the case where 9 and u are arbi- 
trary. The two constants identified in take now the 
form 



sir?9 
1 + cos26l 



[u^ A- ru*) 



E = 



(1 + r^){l + cos^ 



■ u* - 2Lr 



(26) 
(27) 



Like for the Kerr spacetime [4^ a third constant of mo- 
tion V can be identified. For lightlike test particles the 
corresponding first integral is given by 



{u^f (1 -I- cos^9f = V + A cos^9 [L^ - - 



2a / r2 



sin-'e* 



(28) 



For arbitrary 9 the effective potential (to be inserted into 
^ with En = 0) from ([25]) becomes 



V 



off 



2L\l + r^)(- ^^-^ -I- cot^9) 

2{E + 2Lrf V{l + r'^) 



(1 



2(1 -|-cos26')2 



(29) 



In the limit 9 ^ tt/2, V the result (|25|) is recovered, 
to be inserted into ([5]) with En = 2E'^ — 2L^. For time- 
like test particles the first integral receives one additional 
contribution as compared to the lightlike case 



(23) (u^)^(l-Hcos' 



2n\2 



V+Acos^9iL^ 



1 



4L^ 



4M2 sm^9 



(30) 



The constant En on the right hand side of (O is deter- 
mined by En = 2E'^ - 21"^ - 1/{2M^). Circularity of the 
orbits at r = Tc requires that the constants of motion are 
tuned as follows: 



AEL 



l/(2Af2) 



L2 = 



1 



12M2 



(24) 



Therefore, circular orbits are possible only at jrd = oo. 
For L2 = 1/(12M2) the circular orbit is only marginally 
stable because d^V^^ /dr^ = 0. An interesting property 
of timelike geodesies is that they are confined to finite 
radii if 12^2 < but can escape to infinity if 12^2 > 

1/M2. The marginal case leads again to the condition 
L2 = 1/(12M2). 

Lightlike geodesies experience a slightly different effec- 
tive potential. Instead of ([^^ one obtains 



(25) 



= -6L^ - 8LEr 



The constant En on the right hand side of ([5]) is now 
determined by En = 2E'^ - 2L'^. There are no circu- 
lar orbits for lightlike test particles because V^^ = En 
and dV^/dr = cannot hold simultaneously. Lightlike 
particles moving in the plane 9 — tt/2 generically have a 
non-vanishing radial component of the velocity. 



The effective potential for timelike geodesies also contains 
one additional term as compared to the lightlike case 

mi: 



{l + r^)sm'^9 
2Af2(l + cos26i)2 



(31) 



Qualitatively, the potentials (|29p and pTjl behave like 
the simpler potential (|22p : they contain a term quadratic 
in r that can change its sign, a term linear in r whose 
sign is determined by —LE and an r-independent term. 
The discussion about confined geodesies above therefore 
generalizes straightforwardly to generic geodesies. 

We stress that a collection of non-interacting non- 
rotating (L = 0) timelike test particles ("dust") cannot 
be lifted to a regular perfect fluid solution. This can be 
seen directly from the effective potential (j^^ . which be- 
comes positive at sufRciently large r if 12_L2 < 1/Af2. 
Consequently, the radial velocity determined from ([5]) 
eventually becomes imaginary. This observation is just a 
rephrasing of the fact that timelike geodesies with van- 
ishing (or sufficiently small) angular momentum L are 
confined [iTj . In the next sections we circumvent this 
problem by allowing for pressure and find fluid solutions 
that are well-defined everywhere. 
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III. TIMELIKE PERFECT FLUID 

In this section we consider perfect fluid solutions on 
the NHEK background with a timelike velocity 4-vector 



A. Equations of motion 

Without loss of generality we employ the normalization 
u°Ua = -1 (32) 
The perfect fluid energy-momentum tensor 

Tiab a b , D \ ab \ o.b ah , ah ioo\ 

f = pu u F l\ A ■— g +U U (OO) 

must be covariantly conserved 

VaTt = (34) 

Here p is the density, P is the pressure and both are 
connected by the polytropic equation of state 



(35) 



For simplicity we assume for the time being n — I. 

Contracting the conservation equation (|34p with Ub 
leads to the relativistic continuity equation 



(36) 



Contracting the conservation equation ([34)) with leads 
to the relativistic momentum equations 



(p + P) u^Vau" + A"" VaP = 



(37) 



These are our equations of motion. The simple structure 
of the NHEK background (fT5| considerably reduces the 
amount of terms contributing to the equations of motion. 
Moreover, by analogy to the Kerr case [11,13) |4^,[4^ we 
shall exclusively consider velocity profiles that are sta- 
tionary and axisymmetric. 



= M°(r,( 



(38) 



In appendix [A] the equations of motion are displayed ex- 
plicitly for arbitrary velocity profiles of the form above 
((38)) . For vanishing polar velocity — the continuity 
equation ()36p evaluated on the NHEK background ITS)) 
reduces to 

u^'drp +ip + P) drvT = [u*^ = 0] (39) 



equations of motion are still somewhat lengthy we restrict 
our attention to the plane 9 = tt /2 for the time being. In 
additions to the equations of motion the normalization 
condition ((5^ must be fulfilled. The simplest way to 
achieve the latter is by demanding = and therefore 



1 



(40) 



Physically, this velocity profile corresponds to choos- 
ing vanishing angular momentum of the fluid particles, 
see the left Eq. ()23p . The continuity equation holds 
per Ansatz, and the momentum equations evaluated at 
9 = 11/2 lead to only two conditions 



deP = Q 

a (1 -I- r-^) 



[9 = ^/2] (41) 
[9 = ^/2] (42) 



The momentum equations are solved by the following 
density and momentum profiles a,t 9 = ix 12: 



p(r,0 = 7r/2) =po(l + r-')"^ 



P = ap (43) 



Our next goal is to lift the solution in the plane 6 — it/ 2 
given by (|40]) , ([43]) to a solution valid in the whole NHEK 
spacetime. This is indeed possible and our final result is 



1 



u 



M^(l + r2)(l + 



cos^ 



u'' = 

_ _ 



u 



M^(l + r2)(l -hcos20) 

p = po((l + r2)(l+cos20))~^ 
P — a p 



(44a) 

(44b) 
(44c) 
(44d) 

(44e) 
(44f) 



The configuration ((44)) solves the equations of motion 
((36)) . ((37)) . It is gratifying that reasonable values of the 
equation of state parameter, a G [—1,1], always lead to a 
density profile that is finite in the whole spacetime. Note 
also that there are no divergences induced by the depen- 
dence on the polar angle. To the best of our knowledge 
((m is the first exact perfect fluid solution constructed 
on the NHEK background ((T5)) . 

We generalize now the solution ((43)) to the polytropic 
equation of state ((55)) with n 7^ 1 (and n > 0) . By solving 
the momentum equations in the plane 9 — tt/2 we obtain 
a density profile of the form 



B. Exact solutions 

We consider now solutions to the equations of motion 
with vanishing polar velocity = 0. A simple class of 
such solutions can be obtained by assuming = 0, which 
obviously solves the continuity equation ()39p . Since the 



p{r,9 ^7r/2) 



- + Po (1 

a 



(45) 



and P — a p". For negative a the solution is always glob- 
ally well-defined. For positive a and n > 1 the density 
inevitably becomes negative in some region of spacetime 
and such solutions must be discarded. For positive a and 
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< n < 1 the solution is globally well-defined provided 
that the inequality apo > 1 holds. 

The thickness of an accretion disk is usually derived 
from a balance between geodesic deviation (gravitational 
forces) and other physical forces (see for instance the ap- 
pendix of (§1). We calculate now the geodesic deviation 
between two nearby test particles. Given a distance vec- 
tor the acceleration 5° between these particles can be 
calculated from the geodesic deviation equation 



(46) 



where is the 4- velocity. For the solution (1441) and for 
a vector ^ 0, x'^ — otherwise, we obtain 



1 - 3 cos^ 



1VP{1 + cos26')3 



(47) 



Interestingly, there is a zero in the acceleration for a cer- 
tain polar angle, even without the presence of additional 
forces. The critical value is given hy 6 0.955 (about 
54.74 degrees). Curiously, this value of 9 also leads to a 
vanishing Chern-Pontryagin density (|2ip . 

We do not further dwell on the perfect fluid case since 
our main interest are viscous fluids. 



C. Shear viscosity 

In the presence of shear viscosity the energy momen- 
tum tensor for a viscous relativistic fluid receives an ad- 
ditive contribution 



that the momentum equations (|52|) are fulfilled for our 
solution , provided the viscosity function is chosen as 



%/TTr 



[solves momentum equations] (53) 



Here 770 is an arbitrary function of the polar angle 9. 
Technically, the reason for the persistence of the solution 
(|44|) with ([53|) at the level of momentum equations ([52]) 
is that for = = Q the only non-vanishing compo- 
nents of the viscosity tensor are 11'''^ and 11*'" (actually, 
the latter component also vanishes in our case). How- 
ever, the continuity equation (|5ip is violated unless 770 
vanishes. Thus there is no simple way to lift our perfect 
fluid solution ([44]) to a viscous fluid solution. 

It is important to understand why this happens, so 
that a more refined analysis can resolve this problem. As 
we have pointed out, the perfect fiuid continuity equa- 
tion ([55]) holds identically for — = 0. Therefore, 
the contribution from the viscosity tensor to the conti- 
nuity equation would have to vanish by itself. This is 
not the case for the configuration (|44)) with ((53)) unless 
?7o = 0. We conclude that a configuration for a viscous 
fluid should either contain a radial velocity component, 
7^ 0, or a polar velocity component, 7^ 0, or both, 
unless the identity (V(aWf)))(V°M^) = holds. We have 
not succeeded in flnding exact solutions with viscosity 
for the timelike case. However, we are able to construct 
solutions for the lightlike case discussed below. 



IV. VISCOUS NULL FLUID 



(48) 



The viscosity tensor H"** is given by 

H"'' = -2?7 (A^'^A^'* V(,Ud) - - A"'' V^m") (49) 

3 

The function rj parameterizes the local strength of viscos- 
ity. Three important properties of the viscosity tensor are 
symmetry, tracelessness and transversality: 

The tracelessness is a consequence of neglecting bulk vis- 
cosity. The continuity equation generalizing (|36|) is now 
given by 

Va(pu") +PVaM"+n"^VaM6-0 (51) 

The momentum equations generalizing (|37p are now 
given by 

{p + P) u'^Vau'' + A''^ {VaP + V.n^a) = (52) 

In the Newtonian limit only the momentum equations re- 
ceive corrections from viscosity. It may be checked easily 



In this section we consider lightlike viscous fluids 

u'^ua = (54) 

and construct exact solutions. We also add an electro- 
magnetic field. 

A. Physical preliminaries 

Before writing down any equations we want to include 
some input from physics, because this will guide us to 
flnding appropriate solutions of the equations of motion. 
First of all, as explained in the introduction, the lightlike 
case is relevant for physics at or near the ISCO. Second, 
it is a reasonable approximation to consider vanishing 
angular momentum, because frame dragging effects act- 
ing on particles near the ISCO dominate over the relative 
movement on top of the frame dragging effects. From the 
left Eq. (P5)) we may therefore assume that u'^ — —ru*. 
Third, we want to include viscosity and have experienced 
difficulties with such an inclusion if u'" = = (see pre- 
vious section). Therefore, we should refrain from setting 
both of these velocity components to zero. In addition 
to these physically motivated conditions we assume for 
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the time being the polar velocity to be negligible as com- 
pared to the other velocity components. Consequently, 
we employ again the simplification ^ 0. The only 
unknown functions in the 4-velocity are therefore and 
u'' , and they are related by the null condition u'^Ua = 0. 
Putting all these ingredients together leads to the follow- 
ing Ansatz for the 4-velocity: 

u* = u\r,9) u'' = (l + r2)M' it" = u"^ = -rw* 

(55) 

Additional constraints come from our desire to eventu- 
ally match scalar quantities in the NHEK geometry with 
corresponding quantities in the Kerr geometry evaluated 
at/near the ISCO. However, we have seen in section 
III CI that the ISCO at r = M in the Kerr geometry is 
mapped to arbitrary values of the radial coordinate r 
in the NHEK geometry. Thus, a sensible matching is 
possible only if the corresponding scalar quantities are 
independent from r. Since we also require 4> and t inde- 
pendence for our background solution, all scalars in the 
NHEK geometry that allow for a sensible matching must 
be invariant with respect to all four Killing vectors (|16p . 
(fT7|) — they can only depend on the polar angle 9. Thus, 
we assume 



P{0) 



3^ 



(56) 



where for consistency we have imposed the relativistic 
equation of state for a lightlike fluid. 

The energy momentum tensor of a viscous null fluid is 
given by 
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(57) 



The shear viscosity tensor is slightly different as com- 
pared to the timelike case 



H' 



ah 



-2rj (A'^^A^'^ V(,u<i) - ^ A"" V,w^) (58) 



Instead of the properties ([50)1 we obtain 

jjah ^ jjfca jja^ ^ q Uabu'^u'' = (59) 

but not necessarily H°f,u'' = 0. The continuity equation 
in general simplifies to 



u'^daP + UaV.H"'' = 



(60) 



Together with the Ansatz = 0, P — P{0), the continu- 
ity equation (|60p turns into a constraint on the viscosity 
tensor. 



UaVfeH 



ab 







(61) 



The momentum equations can be presented as 

g'''')9aP + 4PVa(M°M^)-hVan°'' = (62) 



(4uV 



Of course, only three of them are independent from the 
continuity equation 



B. Exact solutions 

Our task is now to find an appropriate velocity pro- 
file u^{r,9) and expressions for p{6) and r]{9). The 
derivations are presented in appendix B. Thus, given the 
Ansatz ([55)) . ((56| we can present already the most gen- 
eral exact relativistic viscous fluid solution in the NHEK 
background: 



uo{9) 



u 



(1 + r2)(l -t- cos^i 
Uo{9) 



1 + cos2£ 

u^ = - 

P = Po 
P-Ipo 
V = Vo{0) 



ruo{9) 



(l + r2)(l-hcos26l) 



(63a) 

(63b) 
(63c) 
(63d) 
(63e) 
(63f) 
(63g) 



The functions 770 and uq are not independent from each 
other but related by the differential equation 



2 sin 9 cos 9 



1 



'70 



cos^ 



■ uo 



2 cot 9 + sin 9 cos ( 
1 + cos26' 



(63h) 



An interesting property of our solution is transversality 
with respect to velocity 



H°''Mb = 



(64) 



Thus, even though the viscosity tensor (|58p in general 
does not exhibit the transversality property (|64p . it has 
this property for the most general solution (|63p compat- 
ible with the Ansatz (|55p . ([56]) . Therefore, our viscosity 
tensor has all the standard properties ([50|) like for an 
ordinary fluid in the Landau-Lifshitz frame [4^ . Other 
useful properties of the viscosity tensor and the velocity 
profile ((63)) are the following ones: 



VaH"^ = 

VaU" = 
u'^Vau'' = 



(65a) 
(65b) 
(65c) 

^ Uo = const. (65d) 



Thus, the viscosity tensor and the velocity are 
divergence- free (|65ap . (j65bp . The velocity profile is 
geodetic (|65cp . It has vanishing twist (vorticity) if and 
only if Uo is constant (j65dp . In that case the veloc- 
ity profile coincides with the velocity profile of lightlike 
geodesies with L — in ((26|) - ((29)l . The shear of the ve- 
locity profile is non-vanishing in general, V (a^b) 7^ 0, but 
its norm vanishes, V(QMb) V'-'^u'''' = 0, as a direct conse- 
quence of the properties (j65bp . (I65cp . 
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If the velocity profile is known (or can be guessed on 
physical grounds) the differential equation (|63h[) leads 
to a prediction for the viscosity profile. Of course, not 
all solutions that are allowed mathematically make sense 
physically — the velocity and/or viscosity profiles may 
have singularities at certain angles, or viscosity may fail 
to be non-negative in the whole spacetime. The simplest 
solution (|63p that is physically meaningful in the whole 
NHEK spacetime has constant uq and constant viscosity 
r]Q > 0. According to (|65d|) this is the only solution with 
vanishing vorticity. A solution with vorticity is given by 
Mo C!C cos + In tan^ | and constant viscosity rjo > 0. This 
solution has a velocity profile with a logarithmic singu- 
larity at the poles 9 — 0,tt, but is regular and physically 
acceptable otherwise. 



C. Conserved currents 

We construct now conserved currents by following a 
standard procedure. We start with the equations of mo- 
tion ([62]) and contract them with the Killing vectors (fT6|l . 
P?)) . By virtue of the Killing equation ^ (akb)o.i,± = 
we obtain currents 

^o,i,± ^Pkli,± + ^W'ubkl,^^ + fcg,i^± (66) 
that are conserved 



Va Jo^,± = 



(67) 



We discuss now briefly some properties of these currents 
evaluated on the solutions The current related to 

energy flux is given by 



ja 



P51 



Au^ut -t- n° 



(68) 



Depending on the velocity profile function U[j{9) this 
current may be spacelike or timelike in certain regions 
of spacetime. For large values of |r| the current (pS]) 
is spacelike between the two horizontal asymptotes de- 
picted in Fig.dl essentially the "ergoregion" , and timelike 
elsewhere, regardless of the velocity profile. The current 
related to angular momentum flux is given by 



(69) 



For any velocity profile function uq{9) this current is 
spacelike everywhere, except at the poles = 0, tt where 
it becomes lightlike. This result can be seen easily by 
exploiting the property 



n„fcn^ = o 



if a = 



9,^ (70) 



which implies JlJai oc klkai oc g^ip. The remaining two 
currents have no analogue in the Kerr spacetime, and we 
have not found a simple physical interpretation for them. 

Integrating the current conservation equation (j67p over 
some spacetime volume V and using Gauss' law leads to 
integral identities 



dv 



(71) 



Here dV is the boundary of the spacetime volume, n" 
the outward pointing unit normal and 7 the determinant 
of the induced metric at the boundary. Of particular 
interest are boundaries at constant radius, r = tq, so that 
n° — ^/<f^ {drY and -y^p/j = \/~9 \/W ■ If we consider 
as volume V an interval r e [ro,ro + e], with the other 
coordinates arbitrary, then the integral identity (I7ip in 
the limit of e ^ simplifies to* 



Thus, if the local identity 







(72) 



(73) 



holds the integral identity (jTS]) is implied automatically. 

A particularly simple case is the integral identity ((72|) 
for the angular momentum flux Ji . Dropping the trivial 
t and integrations we obtain 



7r 

j d9 sin 61 (1 + cos^e*) a^HJ, = 



(74) 



With the formula for 11^0 from the appendix (jBleP and 
the solution for m* (|63ap we obtain 



^ 4?7oUo sin^g 

^ (1 + COS20)3 - 



(75) 



Therefore, the stronger condition ([75)1 is fulfilled and the 
integral identity ([74]) holds trivially. The inequality in 
Eq. ([75)1 is true for physical reasons: both the viscosity 
function 779 and the time component of the velocity (and 
hence uq) must be non-negative. 

We postpone applications of the integral identities ([7^ 
and a comparison to analog identities in the Kerr case to 
the conclusions. 



D. Electromagnetic field 

To describe a relativistic viscous plasma we introduce 
the abelian field strength 

Fab - VaAb - VbAa - daAb - d^Aa (76) 

in terms of the 4- vector potential Aa- Then the homoge- 
neous Maxwell equations 



(77) 



hold automatically. The inhomogencous Maxwell equa- 
tions 



(78) 



Note that the determinant of the metric is independent from r. 
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are sourced by the 4-current j° (with some convenient 
normahzation) . The full energy momentum tensor 

rpab rpab i rjiab /'70\ 

receives a contribution from the Maxwell field: 



(80) 



We consider here exclusively solutions where the Maxwell 
energy momentum tensor is conserved by itself. 

VaTS' = (81) 

We have found several solutions with vanishing or non- 
vanishing current. The latter all turned out to be prob- 
lematic in the following sense: there is always a part 
of the spacetime where the current becomes spacelike, 
j°'ja > 0. This need not be a generic feature, but it 
is a feature present in all the solutions we have found. 
Therefore, we consider only solutions with j" = 0. 

A particular solution of this type is given by the gauge 
field 

n 

At^B ln(l + r2) - £; (cos6' + Intan^-) (82a) 



~ 2B arctanr 

Ar = 

Ag^O 



(82b) 
(82c) 
(82d) 



The field strength can be decomposed into electric and 
magnetic parts with respect to a 4- vector n"", which is 
usually assumed to be timelike and normalized. Since we 
do not have some preferred timelike vector field available, 
we choose instead n" = u°'. The associated magnetic field 
is given by Ba — ^ n^^bacd F^"^ with the non- vanishing 
components 

APsine - 2Af2sin20 ^ ^ 

The associated electric field is given by Ea = n'^Fba , with 
the non-vanishing components 



sine 2A'P sin^e 



(84) 



The action invariant — | FabF"-'' for the solution ([5^ is 
given by 

_ 1 „ E^ - B^ 

~ 4 



F^.F'^'' 



(85) 



2M4(l-t-r2) sin^6i 

The instanton invariant eabcd F"-^F'^'^ for the solution ([82 
is given by 

8EB 



Tpab Tpcd 

^abcd ±< t — 



(86) 



M2(l + r2) aii^e 

Both invariants exhibit second order poles at the two 
poles 9 — 0,11. The action invariant is zero if £^ = ±i?, 
the instanton invariant is zero if either E or B vanishes. 

Our discussion of electromagnetic fields is far from 
being exhaustive, but it clearly demonstrates that non- 
vanishing electromagnetic fields can be included straight- 
forwardly. 



V. PERTURBATIONS 

In this section we consider small deviations from sta- 
tionarity and axisymmetry. We decompose all fields into 
a background contribution and a fluctuation. Our main 
assumption is that the fluctuations in the viscosity func- 
tion rj are negligible as compared to fluctuations in den- 
sity or fluctuations in velocity, concurrent with previous 
approaches (cf. e.g. [4^). We also assume that fluctu- 
ations do not change the equation of state, P — p/3. 
Thus, we make the Ansatz 



Pf ^ p + e5p 
Pf ^ P + eSP = 
Vf 



■Pf 



€ — 5p 



(87a) 
(87b) 

(87c) 
(87d) 



Quantities with an index / contain both the background 
solution (denoted without index) and the fluctuation (de- 
noted by 5 and scaled by the small parameter e). The 
fluctuation functions Su"" and 5p depend on all coordi- 
nates t, r, 9, (j). We demand that the total velocity be 
lightlike. 



u'iul = 2eUaSu'' + 0{e^) = 



(88) 



The inner product between the background velocity m° 
and the velocity fluctuation vector Su"" must therefore 
vanish to leading order in e. In this work we shall ex- 
clusively consider flrst order expressions in e, and many 
of the equalities below are valid not exactly, but only up 
to higher order terms. Another assumption concerns the 
transversality property (j64p : we demand that this prop- 
erty is maintained to leading order in the fluctuations. 



u'' sm + m Su'' = 



(89) 



The property ((89)) guarantees that also at linearized or- 
der in the fluctuations our viscosity tensor has all the 
standard properties (|50p . For similar reasons as in sec- 
tion |TV3 we require that scalar quantities like pressure 
exhibit no radial dependence. 



dr6P = 



(90) 



We assume further that the fluctuations do not intro- 
duce angular momentum, 







(91) 



The conditions ([88)) and ([91)) imply that the velocity fluc- 
tuation vector has only two independent components and 
conveniently can be decomposed as 



(92) 



The first contribution to the velocity fluctuation vector 
leads to a rescaling of the velocity. We call these 
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fluctuations "scaling fluctuations". The second contri- 
bution gives the velocity profile a ^-component and thus 
provides a vertical component to the velocity vector in 
cylindrical coordinates. We call these fluctuations "ver- 
tical fluctuations". We shall discuss both contributions 
separately. 

Before such a discussion we address an important tech- 
nical issue. Below we shall encounter various partial dif- 
ferential equations of the form 

u'^9a/ = ^ idt + il + r^)dr-rd^)f = (93) 

where / is some physical quantity and is the back- 
ground velocity profile given by ([55)1 . The general so- 
lution of the differential equation (jM]) consistent with 
periodicity in the azimuthal angle is given by 

oo 

f = Y^ f^it ~ arctanr, 0) cos (n(0 + - In (1 + r^))) 

r!=0 

+ ^ f^{t - arctanr, 9) sin {n{(j) + - In (1 + r^))) 



(94) 



The free functions depend on two arguments: the 

first one is the combination {t — arctan r) and the second 
the polar angle 9. 



A. Vertical fluctuations 

We set Su = and keep only 5u^ in the velocity fluc- 
tuation vector ((92)) . 

The 6'-component of the transversality condition ([89]) 
implies 



(95) 



Therefore, du^ must have the same Fourier decomposi- 
tion as the function / in ([94| . The </)-component of the 
transversality condition (|89|) holds identically. The t- and 
r-components of the transversality condition (j89|) are re- 
dundant with each other and yield another condition. 



(96) 



With the definition (j58|) for the viscosity tensor the con- 
dition ([96]) establishes a first order differential equation 
for the fluctuation Su^. 



, n /2m' (1 -3cos26')cot6'\ , „ 

\ Uo 1 + COS^b' 



Its general solution is given by 
uli9) 



6u' = 



sin 9 (l + cos^6 



■/(i,^,<^) 



(97) 



(98) 



Of course, the integration function f{t,r,(j>) is not ar- 
bitrary, but must have the Fourier decomposition (|94p . 



For any background profile that has non-vanishing uq 
at the poles 9 = 0,7r, the solution (^5]) necessarily di- 
verges there. If 6u^ diverges at some points it should not 
be regarded as a small (first order) fluctuation. There- 
fore, without having to solve any of the equations of mo- 
tion wc are led to a rigidity result: there are no vertical 
fluctuations as long as the background velocity does not 
vanish at the poles 9 = 0,7t. This applies to the cases 
1*0 const, and uq oc cos0-|-lntan^ | discussed in section 
IIVBI Vertical fluctuations with a regular profile for Su^ 
can exist only if uo{9) vanishes at the poles 9 — 0,71. 



B. Scaling fluctuations 

We set 6u^ — and keep only 6u in the velocity fluc- 
tuation vector ([5^ . 

The continuity equation (j60[) simplifies to 



u^'daSP + UaVbSn"'' = 
Transversality ([55)1 implies 



(99) 



(100) 



It can be shown that transversality also implies 
MaVf,(5n°^ — 0. Therefore, the continuity equation fur- 
ther reduces to 



u^daSP = 



(101) 



Together with the requirement (j90|) we obtain the follow- 
ing result from the general solution 



SP = SP{9) (102) 
Similarly, from the condition (jlOOp we obtain 

oo ^ 

5u = ^ (5m^ (i — arctan r, 9) cos {n{<j> + - In (1 + r^))) 

oo ^ 

+ ^ Su'^Jt - arctanr, 9) sin (n(0 -f - In (1 + r^))) 



n=l 



(103) 



The momentum equation (j62p with free index </> sim- 
plifies to 



^asm - 



(104) 



and is fulfilled identically. The momentum equation 
with free index 9 then simplifies to 



deSP = 



(105) 



which implies that SP = const. A constant shift of the 
constant background pressure can be absorbed by a re- 
definition of units and therefore no physical pressure fluc- 
tuations are encountered in scaling fluctuations. The ve- 
locity fluctuations, however, can be non-trivial. 
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The remaining two momentum equations with free in- 
dices t, r are redundant with each other. We consider the 
r equation. 



(106) 



It leads to a second order differential equation in the co- 
ordinate 9. Since we can linearly superpose fluctuations 
it is sufficient to consider one of the Fourier modes. We 
take 

Su = /,';(i-arctanr, 6) cos {n{(j) + i In (1 + r^))) (107) 

for some fixed integer n. The second order differen- 
tial equation descending from the r-momentum equation 
(|106p reads explicitly 



2uq (3- cos26l)cot( 



uo 



1 + cos^9 
2(1 + coi 



4sin^6i 



-)/,'; = (108) 



With rjQ and uq given this equation can be solved. 

We consider first the simplest case, uq = const, (and 
therefore 770 = const, as well). Then the general solution 
for the velocity fluctuation is given by 

00 ^ 
5u^^ cos (n(0 + - In (1 + r^))) 

(fn i^ " arctanr) cosh0„ + f^it ~ arctanr) sinh 9„) 
00 ^ 
+ sin {n{<j, + - In (1 + r^))) 

n=l 

[fn'^{t — arctanr) cosh0„ -I- f^'^it — arctanr) sinh0„) 

(109) 



with 



- (cos6 + Intan^-"! 
2 ^ 2^ 



(110) 



where the functions f^'^, Z^", and are arbitrary. 
At the poles 9 ~ 0,n the Fourier modes diverge because 
of the term lntan^| appearing in (|110[) . We should dis- 
card fluctuations that diverge at some points, because 
clearly the linearized Ansatz is no longer valid for such 
contributions. The only admissible fluctuation therefore 
comes from the zero mode, n = 0, and the regular part 
of the fluctuation du becomes independent of the polar 
and azimuthal angles, 



/o {t — arctan r) 



(111) 



The arbitrary function /q is not fixed by any require- 
ment so far. It is remarkable that there can be no 
oscillations in the angular coordinates. If we require 
that the fluctuations vanish at some initial time t — 
to then fo(to — arctanr) must vanish for all values of 



r, which implies that fo must vanish in the interval 
(io — 7r/2, tQ+TT/2). The discussion above points towards 
a rigidity statement: it seems unlikely that there are any 
non-trivial linearized perturbations for a background ve- 
locity profile with Uq — const. 

We consider now generic functions uq, rjo. In that case 
it is quite difficult to solve the second order differential 
equation (jlOSp . Moreover, it may well be that again all 
Fourier modes have to vanish because of regularity re- 
quirements. We focus therefore on the zero mode n — 0. 



6u — A{9)fo{t — arctanr) 
Then pUg)) simplifies to 



(112) 



A" + A'(2^ + ^ + t£^^5!^) = (113) 

V Uq ?70 1 + C0S26' / ^ ' 

where prime denotes differentiation with respect to 9. 
One solution is given by ^ = const. The other solution 
is obtained by defining y := UqTJqA' , which leads to the 
first order differential equation 



(3 - cos^e) cot 9 

y +y — — = 

1 + cos^t^ 

The general solution for the function A is given by 



(114) 



Ai9) =Ao+A, f — 
J ^0 



1 + cos X 



{x)r]o{x) sinx 



dx 



(115) 



If either the velocity profile or the viscosity has a zero for 
some value of 9 the second term in (|115p has a singularity 
and for consistency we must set Ai = 0. Even if uq and 
770 are non-vanishing and regular throughout the whole 
spacetime we may have to set ^1 = for consistency. An 
example is the simple solution 770 = const., uq = const., 
which leads to an integral in the second term of (I115|) 
that diverges at the poles 9 = 0,7t. Typically, the only 
allowed fluctuations have 7^ 0, Ai = 0. Without loss 
of generality we may set Aq — I. Consistent linearized 
perturbations for the zero mode n — are therefore given 

by 



Su'^ — /o(t — arctanr) 



(116) 



By the same token as before we obtain a rigidity result: 
if (Ju" is supposed to vanish at an initial time 1 = 1^ the 
function /q must vanish in the open interval (to~7r/2, io + 
V2). 

From the discussion above we may conjecture that 
there are no physically relevant flrst order scaling fluctua- 
tions. If this is true then the first non-trivial fluctuations 
can emerge only at second order. 



VI. ASTROPHYSICS AND KERR/CFT 

In this section we summarize and interpret our results, 
and put them into the perspective of the diverse litera- 
ture. 
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In section |T] we considered the black hole 
GRS1915+105, whose spin is nearly maximal. The 
dynamics close to the innermost stable circular orbit 
(ISCO) could therefore be well described by considering 
viscous fluid solutions on an extremal Kerr background 
in the near horizon region. A non-standard feature of 
the fluid is that on the ISCO its velocity vector should 
not be timelike but lightlike, since the ISCO coincides 
with the extremal horizon. In section [ll] we introduced 
the near horizon extremal Kerr (NHEK) geometry 
P^ . The NHEK spacetime is obtained by considering 
a double limit of the Kerr spacetime: near horizon 
approximation and near extremality of the Kerr black 
hole. We reviewed geodesies of timelike and lightlike test 
particles in the NHEK background for fixed polar angle 
9 = 7r/2. We generalized these results to arbitrary polar 
angles 9 and polar velocities and identified a third 
constant of motion T) analogous to the Kerr case 

(4^ . In section [hi] we constructed exact solutions for a 
timelike perfect fluid on the NHEK background, 

but we had difficulties in including viscosity. Since the 
NHEK geometry describes the physics on the ISCO, 
physically relevant results can be obtained if we demand 
that the velocity vector be lightlike instead of timelike. 
Consequently, our next step consisted in finding exact 
solutions for a viscous null fluid in section IIVI We 
started from the conservation equation of the energy 
momentum tensor for a viscous null fluid, presented 
in ([FT]) . The results that we obtained for the velocity 
profile, density, pressure and viscosity are presented in 
(I63p . For any given velocity distribution we predict 
uniquely the viscosity function, up to rescaling (j63h[) . 
We added an electromagnetic field ((75)) . Consequently, 
the energy momentum tensor receives an additional 
contribution given by ([50)1 . For the velocity profile ([S5)) . 
one particular solution for the gauge field that does 
not induce any current takes the form ()82p . In section 
[y] we considered first order perturbations around the 
background solution ([55)1 . Our assumptions ([55)1 - ([01)1 
led us to separate the perturbations into scaling and 
vertical fluctuations ([0^ . We obtained some rigidity 
results from which we concluded that there are no 
physically significant first order perturbations. 

In summary, it is generally not possible to find exact 
solutions to the viscous fluid equations in the presence 
of electromagnetic fields. We were able to succeed due 
to the simplifications implied by the highly symmetric 
NHEK background ([T5)) . 

We address now potential applications in astrophysics. 
A quantity of interest in numerical simulations of accre- 
tion disks is the viscous torque (cf. e.g. [l^, [lH; elec- 
tromagnetic torque [49] does not arise for our solution). 
Namely, some boundary condition should be imposed on 
the viscous torque, and it is debatable what is the ap- 
propriate boundary condition at the ISCO. It may seem 
natural to demand vanishing viscous torque at the ISCO, 
see e.g. [l| for a review, but this is not necessarily a good 
condition (see for instance [soj , where it is also shown 



that viscous torque can change its sign). Viscous torque 
in our case is determined by one component of the viscos- 
ity tensor ()Bip . namely by H^, integrated over a surface 
of constant radius. This result can be derived by full 
analogy to the Kerr case [E^- The integrated conserva- 
tion equation for the angular momentum flux Ji consists 
only of one term, namely the viscous torque integrated 
along the ISCO boundary, see ([7^ . Our result ([75)1 shows 
that for the NHEK background the viscous torque can- 
not change its sign. The result ([75)1 also provides the 
appropriate boundary condition for the viscous torque 
as a function of the velocity profile uq{9). It could be 
interesting to implement (|75p as boundary condition for 
the viscous torque in numerical simulations of viscous 
fluids on (nearly) extremal Kerr backgrounds. 

It is not necessarily straightforward to translate phys- 
ical results obtained in the NHEK geometry to corre- 
sponding results in the Kerr geometry. In this context 
it is worthwhile mentioning that the NHEK coordinates 
(|15p were obtained from the Boycr-Lindquist coordinates 
([T)) by virtue of two coordinate transformations, none of 
which involved the polar angle 9. Thus, the ^-dependence 
of physical quantities like velocity, viscosity or electro- 
magnetic fields is susceptible to a direct translation into 
the Kerr geometry. This is not true for radial depen- 
dences. Since the whole NHEK geometry corresponds to 
the ISCO in the original Kerr geometry it is not even 
clear what a radial dependence in NHEK coordinates 
means in terms of the original geometry. Therefore, we 
should trust only statements that involve polar angles, 
like the result for the viscous torque discussed in the pre- 
vious paragraph or the relation ()63h[l between velocity 
and viscosity. 

In our discussion we have encountered some special 
values for the polar angle 9: the values 9 « 0.82 (and 
9 K, 2.32) separate the "ergo-region" where the Killing 
vector dt becomes spacelike from the "normal region", 
see Fig. [1] The current related to the energy flux Jq ([68)1 
has the same asymptotic behavior as the Killing vector 
dt , and it is lightlike for the two critical angles mentioned 
above. The value 9 « 0.955 leads to vanishing Chern- 
Pontryagin density (dDl and vanishing geodesic deviation 
for a timelike perfect fluid ([17)1 . It could be interesting 
to look for experimental signatures close to these special 
polar angles in real data. 

An interesting effect in the rich physics of accretion 
disks are quasi-periodic oscillations (QPOs). As men- 
tioned in section [V) we have not found any QPOs for 
perturbations around our viscous null fluid solution. This 
is consistent with the data selection procedure that led 
to the spin resuhs for the black hole GRS1915-hl05 [1]: 
the first of the three criteria employed by McClintock et 
al. states explicitly that QPOs must be absent (or very 
weak) . 

Of course, all our results are based upon certain as- 
sumptions. We discuss now briefiy which of them could 
be relaxed and how this would influence our analysis. 
Our viscous null fluid solution ([M)) describes only the 
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physics at the IS CO of an extremal Kerr black hole. If 
we drop both assumptions we are back to viscous flu- 
ids on a general Kerr background, which is too compli- 
cated for analytic studies. It may be possible, however, 
to use perturbation theory, taking our viscous fluid solu- 
tion (l63|) as the leading order result. Perhaps the "near 
NHEK" -geometry constructed recently [s^ can provide 
the suitable next-to-leading order background geometry. 
Another important assumption was that the velocity vec- 
tor is lightlike. We motivated this physically by arguing 
that particles at the ISCO of an extremal Kerr black hole 
can only move with the speed of light, but slightly outside 
the ISCO the fluid should have a timelike velocity vector. 
Thus, for certain perturbative considerations the timelike 
fluid may provide a better starting point. However, we 
have not succeeded in finding exact viscous fiuid solu- 
tions, unless the velocity is lightlike. The derivation of 
our exact solution (|63p required some additional assump- 
tions, all of which were motivated. Some of them could 
be relaxed, and therefore we collect here our assumptions 
for future generalizations: we assumed vanishing angu- 
lar momentum; we assumed that the polar component 
of the velocity vector is small as compared to the other 
components, and hence set = 0; we assumed that all 
scalar quantities depend only on the polar angle 9 (f56|) so 
that a matching to scalar quantities in the Kerr geometry 
is possible. Another interesting generalization would be 
the construction of solutions with a globally causal cur- 
rent. We have found only solutions with electromagnetic 
fields where the current is not causal globally or where 
the current vanishes. Finally, our results on first order 
perturbations and our conclusion that QPOs are absent 
are based on the assumptions (|89p and (PT|) . Dropping ei- 
ther or both of these assumptions might lead to different 
results. 

Finally, we address briefly some (rather superficial, 
but expandable) relations to the conjectured (extremal) 
Kerr/CFT correspondence by Guica et al. [2^]. If the 
conjecture holds then the near-extreme black hole GRS 
1915-f 105 is approximately dual to (the chiral half of) 
a conformal field theory (CFT), with left central charge 
cl ~ 2 X 10''^. The Kerr/CFT correspondence starts 
from the same assumptions as the present work: a dou- 



ble limit is taken that leads to the NHEK geometry (|15p . 
The chirality of the dual CFT implies that all physical 
excitations must be lightlike, which concurs with our as- 
sumption that the velocity profile is lightlike. It may be 
worthwhile to recover some of the results above from the 
CFT perspective, like the prediction of the viscosity func- 
tion from the velocity distribution, the specific boundary 
condition for the viscous torque, the apparent absence of 
QPOs or the appearance of some critical angles 9. 
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APPENDIX A: PERFECT FLUID EQUATIONS 
IN NHEK 



Here we collect some useful formulas for timelike per- 
fect fluids, u°'Ua = —1, in the NHEK geometry. We 
assume throughout the paper that the velocity 4-vector 
components depend only on r and 9. With this assump- 
tion the equations of motion simplify considerably. The 
continuity equation (j36[) reads explicitly 



u'-drp + u^dep +{P + p) 
I {?,cos^9 - 1) cot( 



cos^ 



+ drvT + deu'^^ = (Al) 



The momentum equations p7p read explicitly 



u\u'drP + u'^deP) + iP + p)(^ 



f r 

ru u 



4:(cos^9 + 4cos^9 - 1) 
(1 + r2)(l -|-cos2 9) 



u'u' 



, 2 sin 9 cos 9 
1 + cos2 9 



4 sin^ 9 



(l + r2)(l-|-cos2 6l)2 

u'^dru* + u'^deu'^ = (A2) 



u''iu''drP + u''deP) + 



1 + 



2A//2(H-cos2 
-uV(l-Kr2) 



■9.P+(P + p)((^*)2r(l + r2)S^ 



4sin2 6' 
(1 + cos2 6')2 



cos ^6'-h6cos26>-3 

(1 -t- C0S2 6')2 

2 sin 9 cos 



cos2 9 



+ u'-drU'' + U^dgvJ- = (A3) 
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.Sl„.r 



M2(1 + C0S2 

^^^^ l + r2 ^ ^(1 + cos- 



dgP +{P + p)(- {u'f sin e {—^ 

V 1 + C( 



cos2 9 



,cos''6' + 2cos26i + 9, 



2(l + cos2 6')3 ^ 
A.f „ t A.\ sin 6* COS 6* ^„ « o\ 
0) ^ (1 + cos"^ Oy^ J 



(A4) 



cos^e + ecos^^-a 



(l + r2)(l + cos2 61)2 
r 



cote* 



(cos^e- 2cos26i-3) cot 6* 
1 + cos2 9 



(1 + r2)(l + cos2 I 
I — 



+ 27/ V + u'dru^ + u'>deu'>'\ = (A5) 

1 + cos^ 9 J 



Only four of the five equations above are independent. 



APPENDIX B: VISCOUS NULL FLUID 
EQUATIONS IN NHEK 

The shear viscosity tensor for the velocity distri- 
bution ([55)1 takes the following form: 



n 



te 



ri {I + cos^ 9) {l + r'^)deu^ 



(Bla) 



1 9 

n„ = — M^r] 



sm 



1 + cos2 
1 + cos2 



l + r2 



{r{l+r^)dru' -2u') (Bib) 
(3(l + r2)a^w* + 2ru*) 



n 

He 
He 



M% (1 + cos^ef (u*)2 a^((l + r^) u*) (Blc) 

(Bid) 



r0 



2 

-M^ (1 + cos^ei) (9eu* 

,2 



sm 



1 + cos26l 



M^r]{l + cos^9) dr{{l + r^)u*) 



1 + cos2 



■dr{{l + r'')u') 



(Ble) 

(Blf) 
(Big) 
(Blh) 



The remaining half of the components follows from the 
symmetry, tracelessness and projection properties (j59l) . 
With the results (|B1[) the continuity equation (pT|) sim- 
plifies to 



-r]{2ru^ + {l + r^)drU^) = Q (B2) 



By integrating (jB2p we obtain that has the form 



/(g) 

1 -h r2 



(B3) 



For convenience, we define J {9) = 
(IB3I) takes the form 



^ 1+cos^ e i consequently. 



uo(g) 



(H-r2)(l-f cos2 0) ^^^^ 
The t-component of the momentum equations will then 
reduce to (j63hp . The 9 momentum equation requires the 
density to be constant, p — po (|63p . The r and (j) mo- 
mentum equations are redundant. 



[1] M. A. Abramowicz (2008), 0812.3924. 

[2] J. A. Wheeler, Our universe: the known and unknown 

(1967), public lecture. 
[3] V. Frolov and I. Novikov, Black Hole Physics (Kluwer 

Academic Publishers, 1998). 
[4] S. A. Hughes (2005), hep-ph/0511217. 
[5] J. Casares (2006), astro-ph/0612312. 
[6] B. Aschenbach, Astron. Astrophys. 425, 1075 (2004), 

astro-ph/0406545. 
[7] J. Greiner, J. G. Cuby, and M. J. McCaughrean, Nature 

414, 522 (2001). 
[8] E. T. Harlaftis and J. Greiner, Astron. Astrophys. 414, 

L13 (2004), astro-ph/0312373. 
[9] J. E. McClintock et al., Astrophys. J. 652, 518 (2006), 

astro-ph/0606076. 



[10] J. E. McClintock and R. A. Remillard (2009), 0902.3488. 

[11] K. S. Thorne, Astrophys. J. 191, 507 (1974). 

[12] S. Katos, J. Fukue, and S. Mineshige, Black-Hole Accre- 
tion Disks (Kyoto University Press, 1998). 

[13] J. Frank, A. King, and D. Raine, Accretion Power m 
Astrophysics (Cambridge University Press, 2002). 

[14] N. I. Shakura and R. A. Sunyaev, Astron. Astrophys. 24, 
337 (1973). 

[15] W. Kluzniak and D. Kita (2000), astro-ph/0006266. 
[16] B. Coppi, Europhys. Lett. 82, 19001 (2008), 0802.1694. 
[17] J. M. Bardeen and G. T. Horowitz, Phys. Rev. D60, 

104030 (1999), hep-th/9905099. 
[18] J. Bardeen and R. Wagoner, Astrophys. J. 167, 359 

(1971). 

[19] A. J. Amsel, G. T. Horowitz, D. Marolf, and M. M. 



16 



Roberts (2009), 0906.2367. 
0] A. J. Amsol, G. T. Horowitz, D. Marolf, and M. M. [36 

Roberts (2009), 0906.2376. 
1] O. J. C. Bias, H. S. Reall, and J. E. Santos (2009), [37 
0906.2380. 

2] M. Guica, T. Hartman, W. Song, and A. Strominger [38 

(2008) , 0809.4266. [39 
3] H. Lu, J. Mei, and C. N. Pope, JHEP 04, 054 (2009), [40 

0811.2225. 

4] T. Azeyanagi, N. O gawa, and S. Terashima, JHEP 04, [41 

061 (2009), 0811.4177. 
5] T. Hartman, K. Murata, T. Nishioka, and A. Strominger, [42 
JHEP 04, 019 (2009), 0811.4393. [43^ 
6] D. D. K. Chow, M. Cvetic, H. Lu, and C. N. Pope (2008), 
0812.2918. 

7] T. Azeyanagi, N. Ogawa, and S. Terashima, Phys. Rev. [44 
D79, 106009 (2009), 0812.4883. [45^ 
8] J.- J. Peng and S.-Q. Wu, Phys. Lett. B673, 216 (2009), [46 
0901.0311. 

[29] A. M. Ghczclbash (2009), 0901.1670. [47 
[30] H. Lu, J.-w. Mei, C. N. Pope, and J. F. Vazqucz-Poritz, 

Phys. Lett. B673, 77 (2009), 0901.1677. [48 
[31] G. Compere, K. Murata, and T. Nishioka, JHEP 05, 077 

(2009) , 0902.1001. [49 
[32] K. Hotta, Phys. Rev. D79, 104018 (2009), 0902.3529. 

[33] A. M. Ghezelbash (2009), 0902.4662. [50 
[34] C. Krishnan and S. Kuperstein, Phys. Lett. B677, 326 

(2009), 0903.2169. 
[35] I. Bredberg, T. Hartman, W. Song, and A. Strominger 



D. Marolf, and M. M. Roberts (2009), 
W. Song, and A. Strominger (2009), 



(2009), 0907.3477 

A. J. Amscl, 
0907.5023. 
T. Hartman, 
0908.3909. 

R. Penrose, J. Astrophys. Astron. 20, 233 (1999). 

T. Jacobson and T. P. Sotiriou (2009), 0907.4146. 

S. Chandrasekhar, The Mathematical Theory of Black 

Holes (Oxford, UK: Clarendon, 1992). 

S. M. Carroll, M. C. Johnson, and L. Randall (2009), 

0901.0931. 

B. Carter, Phys. Rev. 174, 1559 (1968). 

I. Novikov and K. S. Thorne, in Black Holes, edited by 
D. C. and D. B.S (Gordon and Breach, New York, 1973), 
p. 343. 

H. Riffcrt and H. Herold, Astrophys. J. 450, 508 (1995). 

H. Riffcrt, Astrophys. J. 529, 119 (2000). 

M. A. Abramowicz, A. Lanza, and M. J. Percival (1996), 

astro-ph/9611101. 

L. Landau and E. Lifshitz, Course of Theoretical Physics, 
vol. Vol. 6; Fluid Mechanics (Pergamon Press, 1960). 
P. Rcbusco, O. M. Umurhan, W. Kluzniak, and O. Regev 
(2009), 0906.0004. 

R. Ruffini and J. R. Wilson, Phys. Rev. D12, 2959 
(1975). 

M. R. Anderson and J. P. S. Lemos, Mon. Not. Roy. 
Astron. Soc. 233, 489 (1988). 



